The role of magnetic polarons in ferromagnetic GdN 
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We report an interplay between magnetism and charge transport in the ferromagnetic semicon- 
ductor GdN, pointing to the formation of magnetic polarons centred on nitrogen vacancies. The 
scenario goes some way to resolving a long-standing disagreement between the measured and pre- 
dicted Curie temperature in GdN. It further constitutes an extension of concepts that relate closely 
to the behaviour of ferromagnetic semiconductors generally, and EuO in particular. 



PACS numbers: 75.50.Pp, 75.30.-m, 75.10.-b 

Intrinsic ferromagnetic semiconductors, which offer the 
freedom to dope into specific conducting channels with- 
out destroying their magnetic behaviour, are of obvious 
technological interest. From a fundamental point of view 
they are again interesting; it is a substantial challenge 
to understand the interplay between their magnetic and 
transport properties. Among the earliest discovered, and 
even now the most studied, is EuO, which shows a mag- 
netoresistance as large as thirteen orders of magnitude 
across the metal-insulator transition at the Curie tem- 
perature (Tc) PUS], the largest in any compound. The 
strong transport-magnetism interplay is further empha- 
sised by the enhanced Tc seen in electron-doped sam- 
ples |4j [5]. Debate continues about whether the ferro- 
magnetic transition is homogeneous El [7] or whether 
it involves magnetic polarons nucleating around magnetic 
impurities [3] EMU]. 

GdN, the prototypical rare-earth nitride compound, 
provides a rich set of comparisons with EuO. Divalent 
Eu and trivalent Gd share the same half-filled 4/ shell, 
with S = 7/2, L = 0, and a net moment of 7 /ie- They 
share also the same rock salt structure, their reported 
T c are both near 70 K ITTTIfTH]. and GdN also shows 
a strong magnetoresistance at Tc P~3j [17]. Theoreti- 
cal treatments of GdN [18-21 reproduce the measured 
electronic band features well, showing agreement as re- 
gards its semiconductor nature p~3j [T3] , the direct band 
gap [22H24] . and features in the conduction-band den- 
sity of states [13] [25] . Even more important for device 
development is that GdN, unlike EuO, has a dispersive 
valence band of delocalised nitrogen states, which raises 
the possibility that it can support both n— and p— type 
conduction. Its potential in realistic spintronics has al- 
ready been demonstrated by its use in a spin filter [26] . 

In contrast the ferromagnetic exchange mechanism is 
still poorly understood. The atomic-like nature of the 4/ 
electrons necessitates indirect exchange, but those levels 



lie too far below the Fermi level to call on the third-order 
perturbation theory applied to EuO [4j [27j [28] . The va- 
lidity of a proposed exchange channel involving the ex- 
cited Gd 4/ 8 level is uncertain [29] ; in particular note 
these levels are similarly far from Fermi level [25] . Small 
energy differences are found among various spin order- 
ings determined within the LDSA+/7 treatment that has 
successfully reproduced the band structure, but they lead 
to an estimated Tc below 25 K [30] , far below the exper- 
imental values. There is even a 1980 report [31], raised 
again recently [32 , that GdN is metamagnetic and metal- 
lic rather than strictly ferromagnetic, though no recent 
data support that phase. 

In view of the Tc discrepancy, there have been several 
attempts to explore free-carrier exchange mechanisms in 
GdN. Metallic versions of the problem are in disagree- 
ment with transport measurements [19] [33] , but there is 
the potential that even in a semiconductor, Ruderman- 
Kittel-Kasuya-Yosida (RKKY) exchange might lead to 
an enhancement peaking at 60 K for electron doping near 
10 20 cm -3 [33] . Such a carrier concentration is easily ob- 
tained in the rare-earth nitrides, which show a propensity 
for nitrogen vacancy formation [14] [I5j [17] . However, Tc 
does not peak; it is found to remain stable near 70 K 
for a very wide range of carrier concentrations [35 , ris- 
ing above that value only for very heavily doped mate- 
rial [15] [35] [36] . 

In this letter we address the nature of the ferromag- 
netic transition in GdN using experimental data from a 
series of epitaxial films with systematically varied doping 
levels. The study is motivated by the existence in the 
literature of data showing a foot at the base of the mag- 
netisation curve [14] [23] [35] , data similar to those taken 
on recent films grown in our laboratory (Figure [TJ [37] . 
In samples with heavy doping the magnetisation nucle- 
ates uniformly near 70 K (sample H of Fig. [I]), which 
we will show is aided by mobile charge carriers, but it 
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FIG. 1: (color online) Temperature dependent magnetisation 
of GdN films with carrier concentrations of < 10 18 cm -3 (sam- 
ple L), 3.3 xlO 20 cm -3 (sample M), and2xl0 21 cm -3 (sample 
H). 



Sample 


Rocking curve 


n (300 K) 




FWHM (°) 


(cm" 3 ) 


M 


2.36 


3.3 x 10 20 


Mi 


2.37 


4.4 x 10 20 


M 2 


2.64 


5.7 x 10 20 



TABLE I: Structure and transport properties of the three 
GdN films shown in Fig. [2] The rocking curve was performed 
on the symmetric (111) GdN x-ray diffraction peak. 

drops to another homogeneous transition near 50 K in 
samples with very low doping (sample L). Samples of in- 
termediate doping (sample M) still show a transition to 
a ferromagnetic state near the 70 K Tq reported so often 
in the literature [T3HT71 138] , but the transition is incom- 
plete until the temperature falls below ~50 K. Reduced 
doping then results in a smaller magnetisation between 
70 and 50 K, rather than the reduced Tq as expected for 
a homogeneous carrier-mediated drift in the exchange in- 
teraction strength. 

It is significant that films grown at temperatures per- 
mitting epitaxial growth lie in the high-conductivity class 
(sample H of Fig. [I]); lower- conductivity films have been 
reported only for ambient-temperature grown polycrys- 
talline films. Thus to explore the situation further we 
focus here on three well-ordered epitaxial films (includ- 
ing sample M) that we have managed to grow, under 
careful temperature control, with carrier concentrations 
between samples L and H in Fig. [T] We have investi- 
gated the magnetic and conducting state of the films, 
including Hall effect measurements to obtain the carrier 
concentration (n) and the mobility (/i = l/(enp) with p 
the resistivity and e the electron's charge) as a function 
of temperature. The full data set shows it is magnetic 
polarons centred on charged nitrogen vacancies that nu- 
cleate ferromagnetism at 70 K. 



The three samples (M, Mi, M2) were grown in a Riber 
molecular beam epitaxy system using NH3 as the nitro- 
gen precursor and a Gd solid effusion cell. The 120 nm 
thick GdN films were deposited onto 100 nm thick (0001) 
A1N layers which were themselves grown on (111) silicon 
substrates. Tuning the growth temperature of the GdN 
films between 650° C and 670° C allowed preparation of 
samples with slightly varied residual doping (see Table |l|. 
The GdN growth was always under N-rich conditions, 
with a beam equivalent pressure of 1.9 x 10~ 5 Torr and 
5x 10 -8 Torr for NH3 and Gd, respectively, which leads to 
a growth rate of 0.12 ± 0.01 fim/h. The epitaxial nature 
of the films was confirmed during growth by reflection 
high-energy electron diffraction. The GdN layers were 
capped with GaN to prevent decomposition in air. X- 
ray diffraction measurements show that the crystallinity 
is similar for each sample, with the GdN fully (111) ori- 
ented and displaying 6-fold symmetry as expected for its 
growth on a (0001) wurtzite surface. All films are fully 
relaxed and no evidence of impurity phases is found. Ta- 
ble [I] also presents the carrier concentration of each film 
at 300 K, measured by the Hall effect. In each case the 
carriers are electrons, as expected for carriers originating 
from nitrogen vacancies. 

The in-plane field-cooled magnetisation of sample M 
under an applied field of 250 Oe is shown in Fig- 
ure ga). The magnetisation shows an upturn below 
^65 K, matching the Curie temperature Tq of 65±1 K 
deduced from the inverse of susceptibility in the param- 
agnetic regime [Fig. [2jc)]. In addition, there is a clear 
shoulder in the magnetisation at about 50 K, coincid- 
ing with the separation of the field cooled and zero field 
cooled curves (not shown). Field dependent magneti- 
sation measurements confirm the vanishing of hysteresis 
above 50 K [Fig. |2fd)]. Films M x and M 2 both show 
similar behaviour, although with less pronounced shoul- 
ders at 50 K. As emphasised in Fig. [I] the shoulder is 
not observed in more conducting samples [TTJ [38] . It is 
worth mentioning that this shoulder is strongly reminis- 
cent of the "double-dome" seen in the magnetisation of 
electron-doped EuO j4l[T2]. 

The temperature-dependent saturation magnetisation 
is plotted in the inset to Fig. [2|a), with the data at low 
temperature fitted by the expected Brillouin function for 
(S = 7/2) [39]. Remarkably this extrapolates to a Curie 
temperature of 50 K, substantially less than the 65 K 
temperature at which magnetic ordering sets in, but co- 
incident with the knee in the magnetisation. The clear 
implication is that the ferromagnetic ordering that sets in 
near 65 K differs from the homogeneous low temperature 
magnetic phase. 

Turning to electron transport, the correspond- 
ing temperature-dependent resistivities are shown in 
Fig. [2fb). Both the room-temperature values and the 
temperature dependencies are in agreement with pre- 
vious reports [131 [14], with a trend that is consistent 
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FIG. 2: (color online) (a) Field-cooled temperature dependent 
magnetisation of sample M. Inset: high- field remanant mag- 
netisation. The dashed curve is a fit of the low temperature 
data to a ferromagnetic Brillouin function, yielding an intrin- 
sic Tc of only 50 K. (b) Temperature dependent resistivity of 
GdN films with slightly varied carrier concentrations, (c) In- 
verse susceptibility versus temperature indicating a param- 
agnetic Curie temperature of 65 K. (d) Coercive field versus 
temperature. 



with the carrier concentrations in Table H The data lie 
between the extremes of resistivity in the literature, in 
which less resistive samples p~6j [17] show a positive tem- 
perature coefficient of resistance (TCR) near room tem- 
perature, while more resistive samples typically show a 
strong negative TCR [141135] . All show a peak near 70 K, 
but note that for the most resistive samples in the liter- 



ature (rsj 10 f^cm at 300 K) the peak shifts downward to 
about 50 K [35]. 

The present three films appear to lie right at the ap- 
parent metal-insulator transition found in the literature. 
Film M2, with the highest electron concentration of the 
three samples, has a very small negative TCR at room 
temperature and has a metallic positive TCR from Tc 
down to 5 K. At lower electron concentrations (M, Mi) 
the negative TCR at 300 K is sequentially stronger and 
the amplitude of the peak near Tc increases. The most 
resistive sample M shows a negative TCR at the lowest 
temperatures indicating that the ferromagnetic phase is 
also semiconducting. Furthermore, a shoulder appears 
in the resistivity below Tc at about 50 K that is not 
seen clearly in the other two films. The appearance of 
shoulders near 50 K in both the resistivity and the mag- 
netisation clearly links the charge carriers and magnetic 
exchange mechanism. This systematic behaviour is rem- 
iniscent of that observed in Eu-rich EuO [2], although 
we note that the resistivity peak is substantially sharper 
in the present GdN samples than in EuO of comparable 
resistivity. 

It is notable that the resistivity of film M bears 
a similarity to the temperature dependence discussed 
many years ago for magnetic metals, with an elbow- 
like anomaly on the low-temperature side of a resistivity 
peak. The anomaly generates a peak in the derivative 
dp/dT that is identified as occurring at Tc [40]. The 
more rounded peak in the resistance at higher tempera- 
tures signals a point at which the spin-spin correlation 
length a satisfies 1/a ~ 2/^, where kp is the Fermi 
wave radius. It is however clear that the present system, 
GdN, cannot be interpreted on that basis on two grounds: 
(i) it is semiconducting without a well-defined Fermi ra- 
dius and (ii) the magnetisation indicates a paramagnetic 
Curie temperature at the resistivity peak, rather than 
the derivative peak. We thus seek an explanation within 
an inhomogeneous model. 

In view of the dominant doping role played by nitro- 
gen vacancies (Vn) it is here that we seek an enhanced 
exchange interaction that may nucleate ferromagnetism 
at 70 K. As a guide we use the Vn energy levels calcu- 
lated by Punya et al. [41 ; each vacancy is expected to 
bind two electrons in a singlet state, with a third elec- 
tron bound only weakly below the conduction band in 
the dilute limit. As the temperature lowers the carriers 
begin to freeze out of the conduction band to occupy the 
third Vn level. This lends a magnetic moment of 1 \±b to 
the vacancies but more importantly provides a large lo- 
cal electron density that can mediate exchange between 
neighbouring Gd ions within the localised cloud, thus 
forming a magnetic polaron [42 . Evidently ferromag- 
netism first nucleates around these sites at a temperature 
of 65-70 K, forming a nonuniform magnetic phase with 
the vanishing coercivity seen in Fig. [2^d) . 

The most lightly doped of the present samples (M) 
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shows an activation that would place the third level 
about 6 meV below the conduction band (see below), 
though more resistive films in the literature would place 
it somewhat deeper [35]. The difference is likely due to 
the high Vn density in these films; note that the car- 
rier density of film M corresponds to one electron per 
100 unit cells, separated by ~ 1.5 nm, very close to the 
~ 1.4 nm radius of the Bohr orbit obtained using an 
effective mass m* = 0.15mo estimated from the band 
structure [22] and e = 4 measured by the subgap reflec- 
tivity [43]. It is significant that this is just the Vn den- 
sity at which both the metal-insulator and the 70-50 K 
ferromagnetic transitions occur. Thus at higher carrier 
concentrations, such as in film M2, the carriers overlap 
and the ferromagnetic phase forms uniformly near 70 K. 
Even higher doping levels can mediate a Tc enhancement 
above 70 K [15] [35] [36] . 

The discussion above emphasises the importance of the 
charge carrier dynamics in determining both the trans- 
port and magnetic properties of GdN. We have thus fol- 
lowed the carrier concentration as a function of the tem- 
perature by measuring the Hall resistance Rh- The car- 
rier concentration, n, is then signalled by the slope at 
fields sufficiently strong to saturate the magnetization, 
and thus also the extraordinary Hall effect. The ordinary 
Hall effect then leads to a high-field slope that measures a 
carrier density that is very nearly independent of temper- 
ature in the most conductive film M2, placing this film on 
the metallic side of the metal-insulator transition. In con- 
trast the most weakly doped film M shows a high-field 
slope that implies the semiconductor-like temperature- 
dependent carrier density plotted in Figure |3|a). The 
carriers are activated above 65 K, with an activation en- 
ergy near 6 meV, confirming the semiconductor nature 
of the paramagnetic phase. The freeze-out is halted at 
70 K, as the gap closes and the activation energy falls 
in the ferromagnetic near- vacancy regions, but resumes 
with a reduced activation energy of 1 meV below the 50 K 
transition to a homogeneous ferromagnetic phase. 

The temperature-dependent mobility plotted in Fig- 
ure [3jb) is quite low in comparison with conventional 
semiconductors, and remarkable for (i) the temperature 
independance above 100 K, and (ii) the rapid rise below 
the Curie temperature. The implication is that electron 
scattering is dominated by the magnetic disorder; in par- 
ticular static impurities would continue to limit the mo- 
bility at low temperature in disagreement with the data, 
he mobility decreases on entry into the inhomogeneous 
magnetic phase below ~ 70 K, and there is a mobility 
minimum (scattering cross section maximum) at 50 K, 
exactly as is expected at the Curie temperature for spin- 
disorder scattering [40, 44 . The resistive anomaly near 
65 K thus has contributions from both spin disorder scat- 
tering and a closing majority-spin gap in the ferromag- 
netic state [22] . 

Now summarising the main results of this experimental 
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FIG. 3: (color online) Temperature dependence of (a) the 
electron concentration n and (b) the mobility \i for sample 
M, obtained from the high- field Hall effect. 



investigation, we have examined the magnetic response 
and charge transport in GdN with carrier concentrations 
across the metal-insulator transition. The data show 
clear evidence of a double ferromagnetic transition, with 
a homogeneous ferromagnetic phase only below 50 K. 
The remarkable interplay between magnetism and charge 
carriers in GdN, and especially its doping dependence, 
signals the presence of magnetic polarons that nucleate 
about nitrogen vacancies. The polarons render a major 
volume fraction ferromagnetic at even such a small va- 
cancy concentration as 1%, and lead to the commonly 
quoted 65-70 K Curie temperature. Stoichiometric GdN 
then has a Curie temperature near 50 K, in closer agree- 
ment with theoretical predictions. This scenario is essen- 
tial for understanding the distinct "double-dome" shape 
of the temperature-dependent resistivity and magnetisa- 
tion. Interestingly, our results have similarity in many 
crucial aspects with the most studied ferromagnetic semi- 
conductor EuO, where closely related concepts, and par- 
ticularly the effects of magnetic polarons, continue to be 
discussed. 

We acknowledge funding from the NZ FRST 
(VICX0808) and the Marsden Fund (08-VUW-030). We 
thank S. Granville and S.V. Chong from IRL for technical 
assistance during magnetic measurements, and T. Maity 



5 



for providing the data for sample L. 



* Electronic address: lben.ru ck@vuw.ac.nz 
[1] Y. Shapira, S. Foner, and T. B. Reed, Phys. Rev. B 8, 
2299 (1973). 

[2] M. R. Oliver, J. O. Dimmock, A. L. McWhorter, and 
T. B. Reed, Phys. Rev. B 5, 1078 (1972). 

[3] T. Penney, M. W. Shafer, and J. B. Torrance, Phys. Rev. 
B 5, 3669 (1972). 

[4] A. Mauger and C. Godart, Phys. Rep. 141, 51 (1986). 

[5] R. Sutarto, S. G. Altendorf, B. Coloru, M. Moretti Sala, 
T. Haupricht, C. F. Chang, Z. Hu, C. SchiiBler- 
Langeheine, N. Hollmann, H. Kierspel, et al., Phys. Rev. 
B 80, 085308 (2009). 

[6] P. Sinjukow and W. Nolting, Phys. Rev. B 68, 125107 
(2003). 

[7] T. Mairoser, A. Schmehl, A. Melville, T. Heeg, 
L. Canella, P. Boni, W. Zander, J. Schubert, D. E. Shai, 

E. J. Monkman, et al., Phys. Rev. Lett. 105, 257206 
(2010). 

[8] J. B. Torrance, M. W. Shafer, and T. R. McGuire, Phys. 

Rev. Lett. 29, 1168 (1972). 
[9] M. Arnold and J. Kroha, Phys. Rev. Lett. 100, 046404 

(2008). 

[10] M. S. Hillery, D. Emin, and N.-L. H. Liu, Phys. Rev. B 
38, 9771 (1988). 

[11] C. S. Snow, S. L. Cooper, D. P. Young, Z. Fisk, A. Com- 
ment, and J.-P. Ansermet, Phys. Rev. B 64, 174412 
(2001). 

[12] P. Liu and J. Tang, Phys. Rev. B 85, 224417 (2012). 
[13] F. Leuenberger, A. Parge, W. Felsch, K. Fauth, and 

M. Hessler, Phys. Rev. B 72, 014427 (2005). 
[14] S. Granville, B. J. Ruck, F. Budde, A. Koo, D. Pringle, 

F. Kuchler, A. Bittar, G. V. M. Williams, and H. J. Tro- 
dahl, Phys. Rev. B 73, 235335 (2006). 

[15] K. Senapati, T. Fix, M. E. Vickers, M. G. Blamire, and 

Z. H. Barber, Phys. Rev. B 83, 014403 (2011). 
[16] M. A. Scarpulla, C. S. Gallinat, W. S. Mack, J. S. Speck, 

and A. C. Gossard, J. Cryst. Growth 311, 1239 (2009). 
[17] B. M. Ludbrook, I. L. Farrell, M. Kuebel, B. J. Ruck, 

A. R. H. Preston, H. J. Trodahl, L. Ranno, R. J. Reeves, 

and S. M. Durbin, J. Appl. Phys. 106, 063910 (2009). 
[18] P. Larson, W. R. L. Lambrecht, A. Chantis, and M. van 

Schilfgaarde, Phys. Rev. B 75, 045114 (2007). 
[19] C.-G. Duan, R. F. Sabiryanov, W. N. Mei, P. A. Dowben, 

S. S. Jaswal, and E. Y. Tsymbal, J. Phys.: Condens. 

Matter 19, 315220 (2007). 
[20] A. N. Chantis, M. van Schilfgaarde, and T. Kotani, Phys. 

Rev. B 76, 165126 (2007). 
[21] C. Mitra and W. R. L. Lambrecht, Phys. Rev. B 78, 

195203 (2008). 



[22] H. J. Trodahl, A. R. H. Preston, J. Zhong, B. J. Ruck, 

N. Strickland, C. Mitra, and W. R. L. Lambrecht, Phys. 

Rev. B 76, 085211 (2007). 
[23] H. Yoshitomi, S. Kitayama, T. Kita, O. Wada, M. Fuji- 

sawa, H. Ohta, and T. Sakurai, Phys. Rev. B 83, 155202 

(2011). 

[24] H. Yoshitomi, S. Kitayama, T. Kita, O. Wada, M. Fuji- 
sawa, H. Ohta, and T. Sakurai, Phys. Status Solidi C 8, 
488 (2011). 

[25] A. R. H. Preston, B. J. Ruck, W. R. L. Lambrecht, 
L. F. J. Piper, J. E. Downes, and K. E. S. H. J. Tro- 
dahl, Appl. Phys. Lett. 96, 032101 (2010). 

[26] K. Senapati, M. G. Blamire, and Z. H. Barber, Nature 
Materials 10, 849 (2011). 

[27] T. Kasuya, I.B.M. J. Res. Develop. 14, 214 (1970). 

[28] V.-C. Lee and L. Liu, Phys. Rev. B 30, 2026 (1984). 

[29] T. Kasuya and D. X. Li, J. Magn. Magn. Mater. 167, LI 
(1997). 

[30] C. Mitra and W. R. L. Lambrecht, Phys. Rev. B 78, 
134421 (2008). 

[31] P. Wachter and E. Kaldis, Solid State Commun. 34, 241 
(1980). 

[32] P. Wachter, Results in Phys. 2, 90 (2012). 

[33] C.-G. Duan, R. F. Sabiryanov, W. N. Mei, P. A. Dowben, 

S. S. Jaswal, and E. Y. Tsymbal, Appl. Phys. Lett. 88, 

182505 (2006). 

[34] A. Sharma and W. Nolting, Phys. Rev. B 81, 125303 
(2010). 

[35] N. O. V. Plank, F. Natali, J. Galipaud, J. Richter, 

M. Simpson, H. J. Trodahl, and B. J. Ruck, Appl. Phys. 

Lett. 98, 112503 (2011). 
[36] K. Senapati, T. Fix, M. E. Vickers, M. G. Blamire, and 

Z. H. Barber, J. Phys.: Condens. Matter 22, 302003 

(2010). 

[37] Sample L was grown using the methodology described 
in Ref. [14] but under a higher N2 partial pressure of 
4 x 10 -4 mbar. Sample M is the present work, and the 
growth conditions of sample H can be found in Ref. [38] . 

[38] F. Natali, N. O. V. Plank, J. Galipaud, B. J. Ruck, H. J. 
Trodahl, F. Semond, S. Sorieul, and L. Hirsch, J. Cryst. 
Growth 312, 3583 (2010). 

[39] S. Blundell, Magnetism in Condensed Matter (Oxford 
University Press, 2001). 

[40] M. E. Fisher and J. S. Langer, Phys. Rev. Lett. 20, 665 
(1968). 

[41] A. Punya, T. Cheiwchanchamnangij , A. Thiess, and 
W. R. L. Lambrecht, Mater. Res. Soc. Symp. Proc. 1290, 
10.1557/opl.2011.383 (2010). 

[42] J. M. D. Coey, M. Venkatesan, and C. B. Fitzgerald, 
Nature Materials 4, 173 (2005). 

[43] M. Azeem, private communication. 

[44] P. de Gennes and J. Friedel, J. Phys. Chem. Solids 4, 71 
(1958). 



